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The general expression for the second order large scale curvature perturbation in the form of a 
functional over a background solution is derived. The explicit expressions was obtained for two 
special forms of the inflationary potential. In the considered cases, it is shown that a significant 
level of non-Gaussianity can be generated during the super-Hubble evolution only if nonadiabatic 
perturbations are non- negligible at the end of inflation. 
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I. INTRODUCTION. 



o 

Primordial non-Gaussianity has emerged as one of important probes of the inflationary era. In single field slow-roll 
inflation, the non-Gaussianity is suppressed by the slow-roll parameters [H, 0], but it can be sizeable in models 
with a transient violation of standard slow-roll (for example, due to a step feature in the inflationary potential (J|). 
[T i i Non-adiabatic perturbations produced during multi-field inflation may lead to generation of detectable deviations from 
■ Gaussian distribution after inflationary stage. For instance, it can take place at an inhomogeneous end of inflation 
£S) ! @> 0> 0' a * * ne reheating Q, [To| or in the curvaton scenario fill ) . 

It has also shown that significant non-Gaussianity can be obtained during slow-roll multi-field inflation [12j , [Hj|. An 
important feature of most of the known examples of this type is that large scale curvature perturbation evolves even at 
the end of the inflationary stage (see, however, [HI])- Hence, such models are incomplete without an understanding of 
the evolution of the cosmological perturbations until they become nearly adiabatic fl5| . In particular, it is necessary 
to consider the non-Gaussianity behavior at reheating. This raises the question about principal possibility of slow-roll 
models with large non-Gaussianity and adiabatic spectrum. 

Most computations of the non-Gaussianity in multi-field inflation have been carried out in the framework of the 
^-formalism [3, [Tzj], 0, El, H3], 

using labeling the trajectories by slow-roll integrals of motion [2l|, [2 
In an attempt to go beyond standard assumptions, some extensions of the <5A/"-formalism are worked out [2 



fa, 

0. Also, the covariant formalism (27|, 0, [2<j, |3| and the long- wavelength formalism (Ml, [13|, [331^01 are 
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developed. Although all approaches mentioned above are geometrically transparent, here we use method of employing 
\ the slow-roll Klein-Gordon equations. 

The structure this paper is as follows. Chapter ITTI describes the model and basic equations. In Chapter HID slow- 
£^ | roll equations are written and general slow-roll expression for the curvature perturbation in two-field inflation is 
. obtained. For cases of product and sum potentials, the corresponding expressions are written explicitly. In Chapter 
" IIV1 using decomposition of perturbations on the adiabatic and entropy ones, it is shown that large values of /jvl can 
\& . be generated only at sufficiently large entropy perturbations at the end of inflationary stage. We summary our results 
i ' m Section El 

O ' 
CO 



II. MODEL AND BASIC EQUATIONS. 

Let us consider two canonical scalar fields ip and x minimally coupled to gravity 



The Lagrangian density is given by 



£ = -\<P*<Pw ~ Ix'^X^ - U(<p, X ). (2.2) 



'Electronic address: koshna71@inbox.ru 



2 



The Hilbert energy-momentum tensor is denned as 



dC 

T" v = -2- — + gT£, (2.3) 



which gives 



Tt = <pP<p. u + r M X;„ - K ( £%, X) + 2^ + ^"Xyx ) ■ (2-4) 



The energy-momentum tensor (|2.4[) . at least up to second order perturbations, can also be treated as energy- 
momentum tensor of a perfect fluid 

T» v = (j> + P)u"u v + P8» v , (2.5) 

where p is density, P is pressure and is 4-speed. 

We split any quantity / into a homogeneous background and small inhomogeneous perturbations 

f{ V ,x*) = fWW+SfWfax') + \6f {2) {Tl,j) + ... • (2.6) 

where r\ is conformal time, the Latin index i takes values from 1 to 3, and indices in brackets indicate the order of 
perturbations. 

Including second order perturbations, the line element around a spatially flat Friedmann-Robertson- Walker back- 
ground has the form [35| 



ds = -a 



fa) {(1 + 2<(>M + <pW)dr) 2 + + B^dridx* + [(1 - 2^ - </> (2) )<% + 2E$ + iV} . (2.7) 

The metrics perturbations can be classified into scalar, vector, and tensor types according to their transformation 
properties under spatial coordinate transformations on constant-time hypersurface [36] ■ Here we will consider only 
the scalar type perturbations, i.e., we set 

B (i) =B W B V=B?, E^=E%, E%=E%. (2.8) 

The freedom of coordinate choice can be used to impose gauge constraints . For example, one can use the uniform 
density gauge (<5pW = Sp^ = 0) or the uniform curvature gauge (tp^ = ip^ = 0). 

In the analysis of perturbations, a special role is played by the quantities that conserved on large scales if the 
pressure perturbation is adiabatic. The first conserved nonlinear gauge-invariant quantity was obtained in ref. [35| . 
The commonly used quantity is curvature perturbation £ on uniform density hypersurfaces, which is ambiguously 
defined in nonlinear case. We adopt the definition of refs. [l|, [37(, [38j. On large scales, this quantity is associated 
with the variable of Malik and Wands I3a| 



A 1 ) -_ ,/,U) 



-MW 

P(1)=P(2)=0 



A 2 ) --nhP* 



(2.9) 

P(1)=P(2)=0 



by relations (37|, [38 



Gauge transformations allow to write the second order curvature perturbation A 2 ) in the gauge-invariant form j3?J 

C<*> = ~V> (2) - ^-h {2) + 2 T ^5pU'5pU + 2^V>' -( U fk- H )(^]\ ( 2 .ii) 
P(o) W») P(0) V p (o) J V ^(0) / 

where % = a' /a and prime denotes the derivative with respect to conformal time 77. Most simply, this expression 
appears in the uniform curvature (UC) gauge 

^^-^a^j^s^M-fn^-u) ( 5 M)\ (2.12) 

P(0) iP(o)) \ P(o) ) \ P(o) ) 

In the following, all calculations done in the uniform curvature gauge with conditions E^> = E^ = 0, what fix a 
residual freedom of gauge transformations. In order to avoid cluttering the expressions with too many subscripts, we 
omit the index " UC" below. 
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A. Density perturbations. 



Let us write only a few components of the two-field energy-momentum tensor in the uniform curvature gauge. 
Without imposing gauge conditions, all components of the energy-momentum tensor up to second-order terms can be 
found in Ref. [39J. 

The perturbation expansion of energy-momentum tensor (|2.4p yield 



rpQ 

J (o)o 



1 

2^2 



/ 2 I , 2 TT 



(0) 



T. 



(0)j 



01 (1)0 
AT 

01 (2)0 



2^(0) +^(0) "^(0) 



</(o) + X( 0) J - V(o) V(i) - X(o) 5 X(i) - aU^Stp^ - a 2 t/, x (5x ( i) 
(/(o) 2 + X(o) 2 ) 0(2) " ^(o) V( 2 ) - x{o)Sx'(2) ~ a 2 U tV 6ip( 2) - a 2 U, x S X{2 ) - V(i) 2 - <*X(i)' 
~ <bC(i),i<bC(l) + (v'(o) 2 + X(o) 2 ) B (i),i B (i) + 4 (V(o) <V(i) + X(o)<*X(i)) 0(i) 



(2.13) 
(2.14) 
(2.15) 



a 2 U tVV Sip 2 {1) - a 2 U, xx 6x 2 {1) ~ 2a 2 U \ vx 5ip {l) 5x(i) - 4 (V (0) 2 + x' (0) 2 J </> 2 



"tf (•P(0)*V(i),i + X / (o)*X(i),i) - 



(1) 



AT 
0J (1)0 



where we use the shorthand U w = 2^21, U y = 2^1, J7 WW = a 3 ^ Q> — , £/ YY = a 8 ^ 0) — , U wv = a 8 — . 
The alternative notation (|2.5p leads to expressions 

T (0)0 = ~P(0), 

T (o)j = P (o) Sl j> 

ST (2)0 = - 5 P(2) - 2 (p(0) + P(0)) «(l),i«fl) ~ 2 (p(0) + f(0)) 

where v is the velocity potential. 

The comparison of equations ([2~T3]) . ([2~bi"jl . ([2~T7|) and (l2~T5|l . ([27TU|1 . (j!T2"2"jl gives 

P(0) = i^ 2 " 1 " 2^ 2X (°) 2 + ?7(0) ' 

D 1 ' 2 , 1 ' 2 TT 

P (0) = ^(0) +^(„) -Efo, 

¥?(0)<V(i) +X( )^X(i) 



«(i) 



= —B 



(i) 



<P(0) 2 + X( ) 2 



From equations (|2~T5)l . (f^TTBj) and (|2~2U)1 . ([2~2T]) . it follows now that 



^(!) = ^2 _ (^'(o) 2 + X(o) 2 ) 0(i) + P(o) V(i) + x{o)Sx{i) + a 2 U^5tp (1) + a 2 U tX S X (i) , 
S P(2) = y?(o)<V(2) + X(o)<*X( 2 ) + a 2 C / ^<5<^(2) + a 2 U, x 5x(2) ~ (V(o) 2 + X( ) 2 ) 0(2) 

- (^'(o) 2 + X(o) 2 ) (^(1)^^(1) - 40 2 1} ) + <V(i) 2 + ^X(i) 2 - 4 (<P( )<V(x) + X(o)^X(i)) 0(i) 
+ a 2 U tVV Stpf^ + a 2 U iXX 6xfij + 2a 2 C/ vx (%i ) 5x(i) - ^( )^(i),i + X(o) 5 X(i),i) 
/(o)X' ( o)<V(Vx(i),* ¥'( 0) 2 -X / (o) 2 fx t i\ 



I 2 , / 2 

P(0) +X (0) 



/ 2 . / 2 

¥> (0) +X( 0) 



(2.16) 
(2.17) 



(2.18) 
(2.19) 
(2.20) 
(2-21) 
(2.22) 



(2.23) 
(2.24) 

(2.25) 
(2.26) 



(2.27) 



B. Perturbed Einstein equations. 

At first order, the — component of the Einstein equations can be written as 

2a 2 U {0) <f> {1) + yjV(i) + x'Sx'i + a 2 U. v S ni) + a 2 U. x S X (i) + ^q AB W = °> ( 2 - 28 ) 

where A is the Laplace operator. 
The — i component is 

U(t> {l) = 4vrG^ (0) ^ (1) + 47rG X ' (0) <5x(i)- (2.29) 
The consideration of the off-diagonal i — j components gives 

B' {l) +2HB {1) +<t> {1) =Q, (2.30) 

From the trace of i — j components we find 

H4/(i) + (2W + H 2 )<f> (1) + ^A ( f )il) +2HAB il) +AB{ 1) 

= AttG (-a 2 U, v 5ip (1) - a 2 U, x 5 X{1) + <p{ 0) <V(i) + X(o)^X(i) - (V(o) 2 + X(o) 2 ) 0(1)) ■ ( 2 -31) 
Using equation (I2.30P and background Einstein equations, the latter can be reduced to 

H(f>' (1) = 4^G (-a 2 U^5 ni) - a 2 U x S X{1) + <p{ 0) <%) + X(o)<*X(i) - 2t/ (O )0(i)) • (2.32) 
At second order, the — i Einstein equation is given by 

Hcj)(2),i ~ 4"H0(i)<^(i),i + 2HB^) M B^ + B(i^ ki cp^ - AB^^ 

= 4irG (p[ 0) 5ip2,i+x{o) S X(2)a + 2S(p[ 1) 6ip {1 y l + 28x{ 1) Sx{i),i • (2.33) 
We then use the first order — i equation (|2.29p and take the trace of (|2.33p , which gives (JlJ 

0(2) + = 20( 1} + (^(0)^(2) + X(o) (5 X(2)) - ^A- 1 (^> {1)M B {1 f l - AB (1) A0 (1) ) 

+ ^T A l (Md A ^(D + ^'(i) A6 *W + <V(i),fe<^(if + $x{ lhk Sx {1 ) k ) , (2-34) 
where A -1 is the inverse Laplacian, A _1 (A/) = /. 

III. SLOW-ROLL. 

As usual, we introduce the two-field slow-roll parameters 

v - 16^G \U W J ' x ~ 16ttG \U {0) J ' Vvv ~ 8ttG U {0) ' %x " 8ttG U {0) ' Vxx ~ 8ttG U {0) ' 1 J 
The slow-roll condition 

max{e v ,e x , |t? w |, \rj vx \, \r] xx \} < 1 (3.2) 

leads to the relations 

<P(0) - W( ) - 0, X('o) - ^X(o) - °> ^^(o) 2 + 2^ x '(°) 2 < U °' ^ 
On large scales, for non-decaying modes of perturbations, there are analogues 

Vex) " - 0, tfxc'i) - ^X(i) - 0. (3.4) 
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V(2) " W(2) - 0> <*X(2) - % - 0. (3.5) 
The basic slow-roll background equations are 

3%cp[ 0) + a 2 U tV = 0, 3Wx(o) + = 0, ( 3 - 6 ) 

^ 2 = ^ 2 C/(o). (3.7) 
Also, the following slow-roll expressions are useful 

P(o)=U (0)l P(o)=-^(^ + ^,x)' ( 3 - 8 ) 
7/ P(o) ? _ a 2 K + 2a 2 U%U, W + 2U, v U, x U iV>x + t/ 2 x £/ xx 

7(0) 3 C/ (0) 3 K + U % ■ { -> 

A. Perturbed slow-roll Klein-Gordon equations. 

At first order, for N scalar fields (j>i with the potential U(4>i). I = 1,...,N, the Klein-Gordon equations can be 
written as [40] 



Stf w + mH' I{1) - A5cf> m + a 2 | V, M j - — y m j | Hk(i) = 0. (3.10) 

Considering the non-decaying modes of perturbations on large scales, these equations are reduced in the slow-roll 

to 

3n5(t)' I{1) + ^2 ( a2 U^ I<j>K - 247rG</4 (o) 0' x(o) ) 8 A K {i) = 0. (3.11) 

K 

For two scalar fields (p and the equations f|3 .6(1 . (|3.7p and (|3.1ip give 

a 2 / J7 i<p C7 )¥ , \ a 2 / f7 V C/ X 



= ^^t" - s m + w. [~u - u <™) 5x ^> (3 - 12) 

SX' (1) = £ - ^) ^ + £ - *tx,. (3-13) 
At second order, it is convenient to go to the Fourier representation 

f{v,x i ) = J d 3 kShe^(ik lX l ). (3.14) 
Then, for two functions / and h the convolution theorem gives 

[/ h] k = J d 3 A/ p /j q (3.15) 

where d 3 X = ^£^S 3 (k i - p i - q*). 

The second order Klein-Gordon equations in closed form have been obtained in ref. [4l| . Here we use their simplified 
slow-roll form from ref. [4l| (with some additional terms) as a starting point: 
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<V/(2)k + 2 ^<V/(2)k + k 2 5tp I{2 )k + Y (a 2 U, VKVl - 247rGV^ (0) v4- (0) ) Si PK(2)k 

K 

= - J rf 3 A |a 2 (u^ lVKVL + ^^-V>^iU <PKtpL ^j Sip K{1)p 5Lp L{1)fl 
{ SttC\ 2 x / 

+ ( ~W I ^ ^/f(0)^A'(l)p a 2 U iVl Y ^K(0)^(l)q + ¥>I(0) 51 a2U ,<PK St PK(l) 
^ n / K \ K K 



K K 



hrG 
H 



I | 2 ^§- 6( P'l(l)p Y <PK(p) S <PK(l) n + 2p 2 <5^/(i) p Y ^(0)^K(l)q 
J { q K K 

piq l +p 2 2 piq l \ (1 q 2 +piq l 



+ ^j(o)£ — q2 ~ ~T~ J 5 ^(i)p^(i)q + [2 W~ J ^(i)p^d)qj I • ( 3 - 16 ) 

These equations can be simplified at small k {k <C %), insofar as in this case it is possible to consider only the 
long-wavelength field perturbations at right hand side. Indeed, the classical field perturbations was generated from 
the quantum field fluctuations at the horizon crossing (42|, [43[, [44|. Hence, at p ^> H, q 3> W, there are negligibly 
small quantum fluctuations instead of the classical perturbations. As a result, (|3.16p is reduced to 



M'( 2 )k + 2 ^M(2)k + Y { a2u ^ K <Pi ~ 247rG^j (0) ^ (0) J <Vif (2 )k 

A' 

+ (~^) ^^(o)^(i)p ^a 2 J7 iV3/ ^^ (0) <5v3 A - (1)q + ^ (0) ^a 2 C/ iip;c ^ K(1)q ^ 

+ ^77— a 2 Y^VK(0) S VK(i)p^ u ,<pK<pi S( PiC(in \ ( 3 - 17 ) 

iC A- J 

One can now go back to the coordinate space. Using equation (|3.5p and the background equations (13.61) . (|3.7p for 
two canonical scalar fields tp, x on large scales we get 



where /1 and /a are 



^(0) , 


) <fy>(2) H 


a 2 
" 3% 


(u, xv - 






J <^X(2) H 


a 2 

" m { 


{u, xv - 


u v u x 



^ + m \U lVV - ^ j t m + ^ - j ^(2) = /„ (3.18) 

«2 / TI 2 \ «2 / TT TT \ 

ty>(3) = fx, (3-19) 



V U (°) ^(0) U (0) J 



U vxx ~ ^U xx - 2U, XV £* + 2E^S- ] £ x 2 1} 1 , (3.20) 
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A ^ -5Z<I". 



u. 



'4) 



2^:XXV ^ ^,¥>X 



c/ 2 £7 \ 

^(0) ^(0) / 



v ^(0) ^(0) ^(0), 



(1) 



(3.21) 



B. Slow-roll perturbed Einstein equations. 



Under the slow-roll condition, the equation (|2.29[) can be rewritten as 



and the equation (|2.31[) takes the form 

4ttG , . , AttG 
-^(o)<V(i) ' 



n 



4^G 



H 



H 



1 U, v , 1 U, x 



(0) 



2 [7, 



(0) 



L7 



fcf(o) + U,<pxX\q) , V(o) + X(o) 



3ft 



3ft<^ 



(0) 



AnG ( 2 U,xxX{ 0) +U, V x<P{ 0) , V / (o) 2 +X , (o ) 2 \ 

-a — + ^«X ( o) ^7J- ) d X(i)- 



3ft 



(0) 



(3.22) 



(3.23) 



It is convenient to write the second-order perturbed equations in the Fourier space. At slow-roll leading order, the 
equation (|2.28[1 gives 



The equation (|2.34[) takes the form 



TJ U 2 + U 2 

1 'a 2 "I '^ X( i )k -3W( 1)k 



6ft t7, 



(0) 



2£7, 



(0) 



(3.24) 



/l U 
d 3 A{pqB( 1)p B (1)q } = --^J^ (2)k 



1 U., 



2 £7, 



(pq)"-A ^(i) P 



u, " 



(0) 

L7 2 + £7 2 



#X(2)k + / d X< 2^( 1)p (/)( 1 ) c 



'(o) 



2u i0) 5 ^ ~ ms ^ J + S I a2 ^u^ 5 ^ - ^ 



'(0) 



'(0) 



1 / g2 _j_ pq ^2 _|_ pq \ I 

+3 ^^Za^ ^ fc 2 Mi) P (5 ^(i)q + k 2 f X(i) P ^X(i)q J >• (3-25) 



C. Solution of the slow-roll Klein-Gordon equations on large scales. 



The system of the slow- roll Klein-Gordon equations can be solved using the approach of [45| . We turn now to the 
new variables 



U, 



u 



(o) 



U- 



(3.26) 



The equations (|3.18[) . (I3.19P can be written as 



x' (y ~ x ) 



U, 



(0) 



a'-llnf^V 



f<pi 



(0) 



£7, 



(3.27) 
(3.28) 



s 



Their consequence 

(x-y)' + 
has a formal solution 



(0) 



(0) 



^/ x ) e -^'>df + 7(2) e'- 



where 7 ( 2 ) is a constant and 



e F = exp • 



in< ^ 



(0) 



X' + In 



[7, 



(0) 



Here To is an arbitrary moment of conformal time. We take To to be the moment of horizon crossing 
for the given mode k. 

The equations (|3.27[) . (|3.28p can be integrated now to give 



V = 



df + a( 2 ), 

dT'+f3 {2) , 



where ot(2)i are constants and 



J(r) 



Hence, the second order longwavelength fields perturbations are 



U, v 



5X(2) 



u v r 


U(o) 


U(0) Jr, 




u x r 




U(0) Jr, 





-fx + In 



U 



<p'e F ^ (J + 7(2)) 



(0) 



(0) 



The consideration of equations p.32[) . (|3.33p and p.30p at t = t* gives 

x* - y* = 7(2)e i?(r * ) = 7(2) = «(2) - /3(2)- 
Similarly, we get the first order longwavelength slow-roll solution 



where OVi), /^m, 7(i) are constants constrained by 

7(i) = "(l) -0(1). 

D. Curvature perturbation in slow-roll. 

At leading order in slow-roll parameters, we have 



(3.29) 
(3.30) 

(3.31) 
t* when k — H 

(3.32) 
(3.33) 

(3.34) 

(3.35) 
(3.36) 

(3.37) 

(3.38) 
(3.39) 

(3.40) 



(3.41) 
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xxX + U tXV tp ) 5\(i) + U tX Sx(i)- (3.42) 

Using the background slow- roll equations, the first order equations (|3.12p . (13.131) and the Einstein equations (|3.22p 
(|3.25p . we obtain a simple expression 

Sp(2)k = U ttp S(f(2)k + U tX Sx(2)k + J d 3 X {U^Stp^pSip^q + U }XX 5x(i) P Sx(i)d + 2U,vx 5( P(i)v S X(i)<i} ( 3 - 43 ) 

>r 

Sp(2) = U^8tp (2 ) + U. x Sx(2) + U. vv 5ip 2 {l) + U lXX Sxfi) + 2U tVX 6<p(t)6x(i)- (3.44) 
The substitution of equations ([3"Hj) . (|3~4^1) . QSHHl into (f2~T2"]) gives 
87rGC/(o) 

u 2 + m 



C(2) = TT2 I 1 ^,^(2) + U X S X (2) 

- U tVV U tX — U :VX U :V U yipip U 2 v + 2U, VX U^U^ X + U ax U 2 x 



V ~ ,x 



rr ,X A~<XX",X "#X",Xrr O >X >Y * >V >A ' ~ iff- ,<p jfj, r 2 

" I XX ~^T ^T^x * (^ + c/, 2 x ) 2 X(1) 

- 2 - -u^r - 2 w^w * ' x ) (3 - 45) 

where formal expression for 5ip( 2 ) and 5^(2) are given by the equations (|3.35p and (I3.36p . 

This equation is the main result of this work. It is follows from (|3.45p , that the first order fields perturbations give 
a local contribution to C(2)i both directly and implicitly through the quantities 6if(2), &X(2)- 



1. Product potential, U{ip,x) = V((p)W(x)- 



The formal expressions (I3.38p . (|3.39p and (|3.35p . (|3.36p for product potential can be integrated to give correspond- 
ingly 



and 



6tp {1 



V 



<><P(i)*, ox(i) 



V( ) v,ip* 



<p<p 



V 

w 







r, V 2 









V i0 )(r) V%; SV V>* ' V [Q) (r) 



\ x 



W 



W 2 
_ ,x 

W 2 



W, x (r) W IQ>,2 , W M (r) W (0> 



2)*) 



(3.46) 

(3.47) 
(3.48) 



Using the background slow-roll equations and equations (I3.46P , (13.47P , (|3.47p , the equation (|3.45[1 can be simplified 
to give 



c (5 



^GV%W 2 0) V(Q) , 
K W (o)+ V (o) W2 x V ^ 



SnGV 2 0) W 2 x W{Qh s 

v 2 wi X{2) * 



8ttGV 2 v W 2 0) 



V 



(ftp* 



^) T '(0)",X 



y2 



8^GU ( 2 0) ^ 2 X 



y, 2 

"(0)* 

w 2 



{V%W 2 + V 2 W\f 



,x> 

2V 2 „wl - w 2 x v (G) v w - vl t r, , i r , 



(0)" y ,x I T72 - ^(D» 



^, 2 m u 



00 n v 2 v w 2 x v^w 2 0) 2v 2 v w 2 x -v vv w 2 x v (0) -w, xx v 2 v w {0) v (0> w i0>x 



(3.49) 



This result is in agreement with ref. (46| . where was considered two scalar fields with non-canonical kinetic terms. 



2. Sum potential, U = V(cp) + W{x). 



At first order we obtain the known result |4E 



C-DW, 



(0) 



C + DV, 



(0) 



where 



C 



V, 



(0)i 



V, 



(0), 



D 



The second order equations (|3.35[) . (|3.36[) give 

C 2 



^(0)c 

v ,9 



2CD 



(o) 



+ 2CD 



V v 2 




U(o) U 2 Q) 


u h) 


( v vv w (0) 








w, xx w 2 x 








( w tXx v m 


w 2 x v (0) 




JJ 2 

U (0) 



d2 (v vv w 2 0) v%w 2 0) w 2 x v 2 ; 



(0) 



u 2 



u 2 



_ W( ) - W^o)* , 

7(2) u (0> 



D 



, , W ,*x V (o) W 2 X V 2 Q) V 2 W 2 { 



u 2 

U (0) 



U 2 

U (0) 



J I 2 



7(2) Y? + ^ 2 )' 



where 



7(2) = 7(2) - L- 



(0)4 



^ 2 



17. 



2CD 



(o)* 



(0)* 



K<p<p*W( ) 4 



W,xx*V(o)* 



u, 



(0)* 



(0)* 



The second order curvature perturbation £( 2 ) is 
8ttG 



C(2) = - 



V 2 + w 2 



v w * {v%{w {Q) + v (0)m ) - w 2 x (v (0) - v (Qh )) - v 2 ^ (y% + w 2 x ) 



V 2 W 2 1-2-^ — ^ :XX ' y 



Mi)* 
V 2 



' V 2 W 2 ( A V, VV W 2 X + W, XX V 2 t ' 

y w * y ~ 2 (v 2 +w 2 x ) 2 

5x1 



+ w, xx * (w 2 x (v (0) + W {Q> ) - V% (W {0) - W {0> )) - w% (v% + w 2 x ) 1 ^ 



ZVlWl 1 - 2 ,x 



+ 8^G V (0 ), - 



XX'.Vit 1 ^(1)* <5 X(1)* 



v 2 + w 2 



V M 



8irG W(o)» - 



V%W [0) -W 2 X V {0) \ S Xi 2)* 



V 2 + w 2 



This result is identical to that found in ref. 1221. 



IV. NON-GAUSSIANITY. 



Let us consider the spectrum P^(k) and bispectrum i?^(ki, k2, ka), defined by 

(Ck 1 Ck 2 >-(27r) 3 < 5( 3 )(k 1 +k 2 )P c (k 1 ), 



(CkXk.Cks) = (2^) 3 5 (3) (ki +k 2 +k 3 ) J B c (k 1) k 2 ,k 3 ). 
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The bispectrum is parametrized in terms of products of the spectrum and the dimensionless nonlinear parameter 

Inl as da 

S c (ki,k 2 ,k3) = ^/ WL (k 1 ,k 2 ,k3)(P c (k 1 )F c (k 2 )+P c (k 2 )P c (k3)+F c (k 3 )P c (k 1 )). (4.3) 
At leading order, we have 

(C kl Ck 2 ) - (((Dk.CfDk.) = (2^ 3 )(k 1 +k 2 )P C(1) (k 1; k 2 ), (4.4) 

(CkiCk 2 Ck 3 ) = 2(C(l)kiC(X)k 2 C(2)k 3 ) + 2(C(l)kiC(2)k 2 C(l)k 3 ) + 7j <C(2)k 1 C(l)k 2 C(l)k 3 ) (4-5) 

Since the variable £k is Gaussian, the Wick's theorem yields 

(C(i)k 1 C(i)k 2 [C( 2 i)]k 3 ) = 2(2nfS^(k 1 +k 2 + k 3 )P C(1) (k 1 )P C(1) (k 2 ). (4.6) 
Combining equations (|4.2[) - (|4.6p . one can write 

6 f ft v vl (C(i)k 1 C(i)k 2 C(2)k 3 ) + c - P- (A 7 , 

7/jVL(Ki,k 2 ,k 3 ) - — -2 ^— , (4.7) 

where c. p. means cyclic permutations. 

In the linear perturbation theory, all perturbations are separated on adiabatic and entropy ones. To do this, one 
can introduce the variables Ssm, San) |48( 

Ssn\ = cosQSxn) — sinQStp/i), San\ = cos@S(pn\ +siaQ8xn), (4-8) 

where 

n ^(0) • n *(o) 
cosB^^-^, smB = ^- L 

a' a' 

and 



Perturbations in form, with <5sm = 0, describe adiabatic field perturbations [48]. Perturbations with Sam = 
represent the entropy perturbations. 

Going to the non-linear case, it is very convenient to use the equality [48| 

C (1 ) = -H*ZL (4.9) 



and to introduce the rescaled variables 



Ssm = (4.10) 



San) 

Sam = H—^p-. (4.11) 

Using these first order variables (we do not employ the nonlinear generalization of ref. (la), the second order 
curvature perturbation can be parameterized as 

C(2) = Sa(2) + CWC(i) _ 2Cga(mSsm + C S sSs 2 m (4-12) 

where C'a&, Cgg-, Cgj are some functions of the background fields, and the quantity Sa2 is zero if the initial fields 
perturbations are Gaussian. Neglecting the initial non-Gaussianity of the scalar fields ip and Xi we obtain 

6. . v) _ r 9r (C(i)k 1 C(i)k 2 [C(i)^(i)]k 3 ) + c. p. (C(i)k 1 C(i ) k 2 [^ 1) ]k 3 )+c. p. 

5 (((1)^(1)^ C(i)k 3 ) + C - P- (CajkiCajkaCfiJk,) + C P- 
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The calculation of the correlation functions yields 



6 



/Ari(ki,k 2 ,k 3 ) = Cs-a + C s % 



if (kQif (k 2 ) + P(k 2 )P(k 3 ) + KjkjKfa) 
P c (k!)P f (k 2 ) + P c (k 2 )P c (k 3 ) + P c (ki)Pc(k 3 : 



P C (k!)P C (k 2 ) + P C (k 2 )P C (k 3 ) + Pc(k 3 )P C (k 1 ) 

where we introduce the notation 

K(k) = (C(i)kS-k). 

Spectra and correlation functions are easy to calculate, using the well-known result [4 



(4.15) 



l)k* 



P* x P* 

=e v k, "X(i)k* = ^ff e xk, 



^fc 3 



^fc 3 



(4.16) 



where e v k and e x k are independent Gaussian random variables with zero mean and unit variance. 

The analysis of the expression (|4.14[) is too complicated for general form of potential U(tp,x), so let us consider 
only some particular cases. 



A. Product potential, U(ip, x) = V(p)W(x)- 



In this case we have 



87tGV (0) W (0) V v W, x fV (ohx W {a) 

- dl P(i)* ~ — ox{i)* 



V 2 W 2 4- V 2 W 2 
V ,V VV (0) + V (0) W ,X 



Sa (1) = - 



8nG 



v 2 wK + v, 2 ,w 2 

y ,ip yv (0) ~ V (Q) vv ,X 



V M 



(4.17) 
(4.18) 



The equation (|3.49p can be rewritten as 

1 



C(2) 



8nG 



V <S>)* 



V 



cos 2 e 



(0), 



W 2 

,x* 



w, 



(0)* 



V 2 



V 2 



2 sin cos 9 



sin 6 

(0)* / 



w 2 



cos 2 e 



w. 



(0)* 



w, x v m ~ 



(i) 



2sin8cose ' xx 
V w (o)* 



W 2 

,x* 



V 



V 2 



Vr 



(0)* 



V 2 



-5a 



(2) 



(4.19) 



where 5a( 2 ) is obtained from <5<7(i) by replacing 5ipn\*, 5x(i)* by Sip( 2 )*, Sx(2)*- 

The equation above shows that the coefficient Caa is of the order of slow-roll parameters, i.e. 



lew I « 1. 



(4.20) 



This inequality and equation (14. 13p indicate that in models with large /jvl the entropy perturbation SSm-\ can not be 
neglected. In addition, it follows now that a significant level of non-Gaussianity |/jvl| > 1 is possible if and only if 



G S cr 



P c (kx)P(k 2 ) + P c (k 2 ) JT(ki) + P c (k 2 )P(k 3 ) + P c (k 3 )P(k 2 ) + P c (kx)P(k 3 ) + P c (k 3 )P(kx 

P C (k!)P C (k 2 ) + P C (k 2 )P C (k 3 ) + P C (k 3 )P C (kx) 



e.55 



P(k!)P(k 2 ) + P(k 2 )P(k 3 ) + P(kx)P(k 3 



P c (k!)P c (k 2 ) + P c (k 2 )P c (k 3 ) + P c (ki)P f (kg) 



> 1. 



(4.21) 
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Consider the case that can be easily implemented 

max{|CW|, |C SS |} < 1. 
Then the mandatory condition of the presence of large non-Gaussianity is 

\K(k)\ > P c (k). 



(4.22) 



(4.23) 



For this class of models, large /jyz, is achieved when two conditions are met: 1) the variables Saru and Ssn\ are 
highly correlated; 2) entropy perturbations dominate over adiabatic ones, -Pi(k) ^S> Ps&(k). 
Quantitatively, the equations (|4.16p . (|4.17p and (|4.18p lead to the expressions 



K{k) = (8ttG) 3 sin9 cos 6 
F c (k) = (8ttG) 3 (sin 2 6 



W :X W, 



(0), 



(0) 



w, x w t 



cos 2 e 



V v V(o)< 
F( ) V v * 



<6{V%W 2 + V 2 W 2 x )k 3 ' 
) 6(^2 W 2 + V 2 W* x )k 3 ' 



The inequality (|4.23p can be rewritten now as 



( w ,x W (°h 



V, v V(o)* 
V( ) V <(p * 



» I tan 61 



This condition can be satisfied for some inflationary models [1 



\W iQ) W, 



icotei 



V<p V(o)* 
V(o) V t<ptf 



(4.24) 
(4.25) 



(4.26) 



B. Sum potential, U = V{tp) + W(x). 
In terms of initial fields perturbations, the quantities <5o"(i) and <5s(i) are 
r (F(o) + W( Q) )V V W, X fScp^ Sx(i)* 



(5S(i) = 8nG 
6<J(i) = 8ttG 



V 2 



w 2 

■x 



V 2 



w 



8irG 



V(o)*<V(i)* W(o)*8X(i)» 



The equation f)3 . 55[) yields 



(4.27) 
(4.28) 



C(2) = 



V 2 +W 2 
8nGU 2 Qh 



V% (W {0) + V i0> ) - W 2 X ( V( 0) - V {0> ) 
8nGUf Q> (W 2 X + V 2 V ) 



V 



W 2 X (V (0 ) + W {0> ) - V% (W (0) - W {0) ,) 



8nGU 2 0> (W 2 X + V 2 V ] 



W 



XX* 



r 2 



(i) 



v,w* {y%(w {0) +v {0> )-w 2 x {v {0) -v {0> ))-v 2 ^ (v 2 v +w 2 x ) 

^GU 2 0> U 2 Q) V 2 W 2 X 

w xx * (w 2 x (v {0) +w {0h )-v 2 v (w w -w (0h ))-w 2 x * (v%+w 2 x 

8nGU 2 Q> U 2 Q) V 2 v W 2 x 



Wt 



(0), 



V, 



(0), 



w 2 x v {0) -v%w {0) 
(wi+v 2 v ) 

w 2 x v (0) -v^w {0) " 
(w? x +v,l) , 



C(i)<%) 



V.UWjo) + V(o).) - w 2 x (v (0) - V (0) .) 
^GU 2 0> U 2 0) V 2 V W 2 X (W 2 X + V 2 ) 
W 2 X (V (0) + W {Qh ) - v 2 v (w (a) - W (0) «) 
^GU 2 0h U 2 0) V 2 v W 2 x (W 2 x + V 2 v ) 

(w {0> (w 2 x + v 2 v ) + w 2 x v (i 



w (0> (W 2 



v 2 



w 



v% + w 2 x 

8nGU 2 0) 



^GU 2 Q> U 2 Q) V 2 W 2 X 

v w w 2 x + w, xx v 2 v 

W 2 X ) 



xx* 

2 

-V; 



^(o)-^(o)) 

. (V(o). (w 2 x + ^ ) - w 2 x v (0) + y 2 ^oj) 1 



8^ 2 0) ^ 2 0) V 2 W^ X 



8ttGU {0) (V 2 



Ss 2 ^ — 5a 



(2)' 



(4.29) 
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This equation leads to the same large Jnl condition (|4.2ip . as in the case of product potentials. For the case of 
sum potentials, we have 



n(k) 



(87rG) 3 U {Q) U i0> V v W tX VJ^ + Wf xt 
V 2 + w 2 x 



'w 2 x v m -v%w (0) 



(8^G) 3 C/ (0) * 



V 2 



- 2 {W 2 X V (0) - V%W m ) (V ioh W t 
The inequality (|4.23D gives 



,2 W 2 



V 2 



V 2 



X* 



V 2 



w 2 



(W 2 X + Vl) . (4.31) 



(4.30) 



w 2 x v (0) 



v%w {0) v.w t 



2 

x* 



wi 



V 2 



V 2 



wi. 



w 2 x v {0} 



> 



w 2 

yy ,x 



V 2 
y ,v 



u (0) v v w a 



'w 2 x v {0) 



vlw {0) v {a> w 2 xr -w (0> v 2 



- V 2 / 
V 2 w 2 



W 2 + V 2 



V 2 

)¥>* 



w 2 



w 2 V 2 

vv (0)* v ,ip* 



V 2 

y ,<p 



wi 



V 2 



w 2 



V 2 



w 2 

,x* 



(4.32) 



Some examples of slow- roll models satisfying this condition are presented in ref. [13j ■ 

In general, the coefficients Cgg-, Css do not have to be small. For instance, if the field \ IS subdominant during slow 
roll and 



W (0) < W(o)* 



< 1, 



W 



V 2 



x v {0) < w {0> 



(4.33) 



then 



Css 



v v w (0> 



V 



<p<p* 



V 2 

y ,<p* 



W M U {0) y87rGU (0> 8nGU i 



2 

(0)* 



V% W(o)*^(o)* 



W 2 

,x 



JJ 2 

U (0) 



8-kGU, 



(0)* 



W 



XX* 



8-kGU [ 



(0)* 



(4.34) 
(4.35) 



The coefficients Cjg-, Css can be valuable due the enhancing factor Vtp/W jX . In this regime, the curvature pertur- 
bation £( 2 ) remains constant and entropy perturbation may be small, P# s (k) -C P&cr(\t). An example of such solution 
in the "adiabatic limit" was presented in ref. [14]. Although the entropy perturbation Sim is decaying, it can not be 
set equal to zero, since then the parameter Jnl gets a wrong small value. 



V. CONCLUSION. 



We presented a derivation of general expression for the second order curvature perturbation £( 2 ) in the form of a 
functional over a background solution. For two special cases, the explicit expressions was obtained that reproduce 
known results. It is shown that, in contrast to the linear perturbation theory, the curvature perturbation depend on 
both the adiabatic and nonadiabatic perturbations. Moreover, a significant level of non-Gaussianity can be generated 
during the large scale evolution only if nonadiabatic perturbations are non-negligible at the end of inflation. Hence, 
to compare inflationary theory predictions at large /jv l with observational data, one need to investigate the evolution 
of the large scale curvature perturbation on post-inflationary stages. 
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